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Abstract 

We examine various properties of the continued fraction expansions of matrix eigenvector 
^ slopes of matrices from the SL(2, Z) group. We calculate the average period length, maximum 

period length, average period sum, maximum period sum and the distributions of Is 2s and 
3s in the periods versus the radius of the Ball within which the matrices are located. We 
also prove that the periods of continued fraction expansions from the real irrational roots of 
x"^ + px + q = are always palindromes. 
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I — I 1 Introduction 

Any real number x can be expressed as a continued fraction of this form: 
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. ^ where uq £ Z and Oj G Z+ for all i > 1. For simplicity, we write x = [oq, ai, • • • , flm " " " ]• We 

^ will be looking at quadratic irrationalities, or real roots of the equation ax"^ + bx + c = with 

5^ integer a,b,c . When a root is irrational, its continued fraction expansion is infinite, and always 

periodic past a certain index. Any rational has a finite continued fraction expansion. 

In the Continued Fractions brochure by V.I. Arnold [4] he proposes this problem: "Let us 

examine matrices ^ ^ , with integral a,b,c,d and the determinant equal to 1. 1. We choose 

from these the ones that define a hyperbolic rotation. There is a finite number of matrices 
whose coefficients are not too large, ie. + 6^ + + < A^. For each such matrix there 
exists an eigenvector y = Xx, for which A is a quadratic irrationality and therefor the continued 
fraction expansion of A is periodic. Take this period and calculate how many ones, twos, threes, 

etc. there are in its period, and then average this for all matrices ^ ^ . That is to say, 

take the number of ones in each period and divide by the number of elements in the period for 
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each matrix. Hypothesis: this relation wih approach a Gaussian distribution as N approaches 
infinity." (paraphrased into Enghsh from Russian. See also [3]). M.Avdeeva and B.Bykovski 
solved another problem posed in this broshure regarding the Gaussian distribution of period 
elements of continued fractions [5] . 

Ik I 

We have written a program that generates matrices M G SL(2,Z) , such that M = 

\m n 

is within a Ball or radius r around the origin for all r < 1000. That is, it generates matrices with 
the coefficients k, I, m, n such that k"^ -\- P -\- m? + n? <r'^ and kn — Im = 1 . We then calculate 
the continued fraction expansions of the matrix eigenvector slopes and find some statistics about 
their periods. Because any quadratic irrationality can be expressed as the slope of an eigenvector 
of such a matrix, all possible quadratic irrationalities are included in our statistics as r approaches 
infinity. 

I'd like to extend a special thank you to Oleg Karpenkov for valuable help with the revision 
process. 



2 Methods 



Our program generates all matrices M S SL(2,Z) in a ball of radius r around the origin, such 
' k l\ 

and A; > for all r < 1000 using the algorithm desribed in 2.1. We find one 



that M 



eigenvector 



m n 
1' 

yi 



for each matrix M using yi 



21 • yi IS real, we nnd the 



period of its continued fraction expansion using a variation of the Euclidean Algorithm described 
in 2.2. If yi is rational, we say that its period is 0. We do not make any calculations for the 

^ ' of M because by Lemma 1, ?/2 has the same coefficients in its period as 



second eigenvector 



2/2 



yi, but in reverse order. We count each matrix with coefficient A; > 1 twice, because it has the 

-k -I' 



same eigenvectors as the matrix M 



-m 



-n 



which we had not accounted for. 



2.1 Matrix Generating Algorithm 

We generate all matrices M G SL(2,Z) in a ball of radius r around the origin such that M 
k ll 

and A; > 0. We do this by generating all k between and r. Then for each k we 



generate all m between — \/ — k"^ and Vr"^ — k"^. For each k,m pair, we find a solution to the 

r /j. I 

Diaphontane equation kn — Im = 1 for l,n. This way we obtain one matrix M 



each k, m pair. To obtain the rest of the matrices, we find Ma 



k I + ak 
m n + am 

and a = —1, —2, —3... such that /c^ + (/ + ak)'^ + + (n + am)^ < r^. 



for 

m n 
for a = 0,1,2,... 
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2.2 Computing the Continued Fractions of Quadratic Irrationalities 

We wish to compute the continued fraction expansion of xq = "o+^ov^ ^ ^ 

using the Euchdean Algorithm. Let oq be the integer part of xq and let x'q be the fractional 
part of Xq. Then let xi = ^ and ai be the integer part of xi. We continue this process finding 

o 

xo,xi,--- ,Xn until Xn = Xi for some i < n. Then the continued fraction expansion of x is 
[ao, ai, • ■ ■ ) o.n] and the period of the continued fraction expansion is [ai, aj+i, • • • , a^-i]. 

Notice that xi = "i+^iV^ for ui,vi,zi £ Z that are easily expressed in terms of uo,vo,zo 
and oq. V. I. Arnold conjectured that if we begin with xq that is a real root oi x'^ + px -\- q 
for p,q £ Z, then the Ui and Zi obtained through this algorithm are always divisible by fj. In 
other words, we are always able to cancel out Vi in this special case. This conjecture held in 
all of our calculations. Furthermore, we found that for the general case where xq is a real root 
of ax^ + bx + c for a, 6, c G Z we are also able to cancel out large numbers from the numerator 
and denominator by finding their greatest common divisors. We noticed that after canceling, 
the coefficients Ui,Vi,Wi, and Zi grow very little throughout the calculation of any continued 
fraction expansion. We compared our algorithm with the same algorithm that does not cancel 
out the common divisors and found that our algorithm has a significantly lower big O running 
time. Without canceling the coefficients grow exponentially throughout the calculation, while 
with canceling they are bounded. 

We also discovered a better algorithm that we have not yet implemented. The idea is to find 
at every i-th. step the equation Aix'^ + BiX + Ci that has the root Xi described above. Then if 
we let Gi be the integer part of Xi, the equation for Xj+i has coefficients ^j+i = Aiaf + Biai + Ci 
, -Bj+i = Bi + 2Aiai and Cj+i = Ai . Hinchin [61 proved that the coefficients Ai,Bi and Ci are 
bounded with the bounds Ci = Ai-i < 2\Aoxo\ + |j4o| + |i?o| and Bi can be defined in terms of 
A, and Ci by Bf = B^ - AAqCq + AAid. 

3 Results 

3.1 Average Period Length vs. Radius 

We computed the average period length of the continued fraction expansions of the slopes of 
matrix eigenvectors for matrices within a ball of radius r around the origin for r < 1000. We did 
this by summing all the period lengths for matrices within a given radius and dividing by the 
total number of matrices with real eigenvectors and real eigenvector slopes within that radius. 
We plotted the average period length versus the radius as shown in Fig. 1. We found that the 
average period length grows as Igr and can be very precisely approximated by —6.11 + 1.92 Igr 
for r < 1000. 

3.2 Maximum Period Length vs. Radius 

We computed the maximum period length of the continued fraction expansions of the slopes 
of matrix eigenvectors for matrices within a ball of radius r around the origin for r < 1000. 
We plotted the maximum period length versus the radius as shown in Fig. 2. We found that 
the maximum period length appears to grow as Igr. For example, for a radius of r = 100, the 
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Figure 1: Average period length of continued fractions associated with matrix eigenvectors for 
matrices within a ball of radius r for r < 1000. It can be very precisely approximated by 
-6.11 + 1.92 Igr. 



expansion of whose eigenvector has the period [1, 1, 1, 1, 1, 1, 1, 2]. 

3.3 Average Sum of Period Elements vs. Radius 

We computed the average sum of period elements of the continued fraction expansions of the 
slopes of matrix eigenvectors for matrices within a ball of radius r around the origin for r < 1000. 
We did this by summing all the elements in the periods for matrices within a given radius and 
dividing by the total number of matrices with real eigenvectors and real eigenvector slopes within 
that radius. We plotted the average period sum versus the radius as shown in Fig. 3. We found 
that the average period sum grows faster than Igr. 

3.4 Average of the Quantity (Period Sum / Period Length) vs. Radius 

We computed the average quantity (sum of period elements / period length) of the periods of the 
continued fraction expansions of the slopes of matrix eigenvectors for matrices within a ball of 
radius r around the origin for r < 1000. We did this by summing the quantities Qi for matrices 
within a given radius and dividing by the total number of matrices with real eigenvectors and 
real eigenvector slopes within that radius. Here Qi is the sum of the period elements divided by 
the length of the period for the i-th continued fraction. We plotted the average of the Qi versus 
the radius as shown in Fig. 4. We found that the average of the quantity (period sum/period 
length) grows as Igr and can be very precisely approximated by —2.70 + 3.71 Igr for r < 1000. 



maximum period length is 8. It occurs for the matrix M 
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Figure 2: Maximum period length of continued fractions associated with matrix eigenvectors for 
matrices within a ball of radius r. The maximum period length grows as Igr. 
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Figure 3: Average sum of period elements of continued fractions associated with matrix eigen- 
vectors for matrices within a ball of radius r. The average period sum grows faster than Igr. 
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Figure 4: Average of (period sum/period length) for continued fractions associated with matrix 
eigenvectors for matrices within a bah of radius r. It can be very precisely approximated by 
-2.70 + 3.71 Ig r for r < 1000 



3.5 Maximum Sum of Period Elements vs. Radius 



We computed the maximum sum of the period elements of the continued fraction expansions 
of the slopes of matrix eigenvectors for matrices within a ball of radius r around the origin for 
r < 1000. We found that for a radius r > 5 the maximum sum of the period elements is r — 2 and 



occurs for the matrix M 

l-r+-y/(l-r)2-4 



1 

-1 1 -r 



The eigenvector ei 



of M has the slope Ai 



[l,r-3]. 



which has the continued fraction expansion [—1, l,r — 3, • • •] with the period 



3.6 Appearance of Is 2s and 3s vs. Radius 

The Gauss-Kuzmin distribution gives the probability distribution of the occurrence of a given in- 
teger in the periods of the continued fraction expansions of arbitrary real numbers. . Bykovski 
and Avdeeva proved that this is also true for arbitrary quadratic irrationalities |5l. The per- 
centage of Is, 2s and 3s in the periods of the fractions that we calculated should follow the 
Gauss-Kuzhmin distribution 



Pr(i^ = fe) = -log2 1 



1 



So we would expect that the percentage of Is would tend to 0.415, the percentage of 2s would 
tend to 0.169 and the percentage of 3s would tend to 0.093. It is evident from Fig. 5, 6 and 7 
that the radius is too small for us to see this distribution. 
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Figure 5: Percentage of Is in the periods of continued fractions associated with matrix eigen- 
vectors for matrices within a bah of radius r for r < 1000. The Gauss-Kuzhmin distribution 
requires that the percentage should tend to 0.415 as r tends to infinity. 
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Figure 6: Percentage of 2s in the periods of continued fractions associated with matrix eigenvec- 
tors for matrices within a ball of radius r for r < 1000. The Gauss-Kuzhmin distribution requires 
that the percentage should tend to 0.169 as r tends to infinity. 
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Figure 7: Percentage of 3s in the periods of continued fractions associated with matrix eigenvec- 
tors for matrices within a bah of radius r for r < 1000. The Gauss-Kuzhmin distribution requires 
that the percentage should tend to 0.093 as r tends to infinity. 



4 Palindromes Proof 

Theorem 1. Let xi be a real root of + px + q for p, g G Z. And let [ag, oi, • • • , On— ii o-n 

1 he 

the period of the continued fraction of xi. Then the period is a palindrome in the sense that if 
we look at the period as a cycle, the cycle is symmetric about some point (which might be an 
element of the cycle or be between two elements of the cycle). This conjecture was stated by 
V. I. Arnold fM/ and proven by Francesca Aicardi fJ!^. Here is an alternate proof that we found 
independently. 

Proof. This proof rests on two big ideas. First, if xi and X2 are real roots of 

x^+px + q (1) 

for p,q G Z, they both have the same period. This was proven by V. I. Arnold [2j. Second, 
Proposition 1 states that if yi and 2/2 are real roots of any quadratic equation with integer 
coefficitents , they have the same coefficients in their period, but in reverse order. Because xi 
and X2 are roots of (1), they must have the same period and the period of xi must be the reverse 
of the period of X2. Consequently, the period of xi must be a palindrome. □ 

Proposition 1. Let yi and 2/2 be real roots of ay^ + by + c for a,b, c ^ Z . If we look at the 
period as a cycle, then the period ofyi has the same coefficients as the period of y2, but in reverse 
order. 



Proof. Lemma 1 states that if we let M E SL(2, Z) be a matrix with eigenvectors 
then the period of ei, [ai,a2) •••jfln] is the reverse of the period of 62, [a„, ...02,01]. Lemma 2 



"l" 




"1" 


and 






.^2_ 
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states that we can construct a matrix M G SL(2, Z) such that 
of M. It follows that the period of yi is the reverse of the period of y2- 



"1" 




"1" 




and 









are the eigenvectors 



Lemma 1. Let M be a hyperbolic matrix in SL(2,Z) with eigenvectors 
the period of e\, [ai, 02, an] is the reverse of the period of 62, [a^, ...02, ai]. 





and 


"1" 






.^1. 







□ 

Then 



Proof. The proof is based on the geometric interpretation of a continued fraction as the boundary 
of the convex hulls of integral points in the angles formed by two lines. The convex hull formed 
between the two eigenvectors of a matrix is invariant under a linear transformation of that 
matrix. The boundary is shifted by the matrix along itself. The region of the period on the 
boundary is shifted to another region of the period on the boundary. 

Let Ai and A2 be eigenvalues of M. Because M is hyperbolic and has a determinant of 1, 
we can assume WLOG that |Ai| > 1 and IA2I < 1. . Therefore, by applying M finitely many 
times, we can shift the period of ei onto the period of 62. Consequently, their periods must be 
the reverse of one-another. 

□ 

Lemma 2. Let yi and 2/2 be real irrational roots of ay'^ + by + c for b, c E 7.. Then we can 



construct a matrix M G SL(2, Z) such that 



"1" 




"1" 


and 






.2/2, 



are the eigenvectors of M. 



Proof. The continued fraction expansion of yi has some period. It is known that there exists a 
matrix M G SL(2,Z) which shifts one period of yi into the next one and ' ^ 



yi 



is an eigenvector 



of M. We wish to shows that 



is the other eigenvector of M. Let Ai and A2 be eigenvalues 



of M. Let u = a + fi^ where 7 is the discriminant of the characteristic equation of M and a 
and /3, are rational. We define conjugation for -u as u = a — /?y^. Then y\ satisfies [M — \\L\ 



1 

y\ 

thus 



=0. Consequently it is also true that [M — \\L\ 
is an eigenvector of M, as we wanted. 



0. Thus [M - A2I] 



0, and 



□ 



Remark 1. Ln all of our calculations the greater root ofx^+px + q the period was always either 

exactly a palindrome or would be a palindrome with its first or last element removed. However, 
this is not always true for the periods of the lesser root. We also noticed that for \p\, \q\ < 100 
the period of the first root is never shifted by more than 3 places from that of the second root. 
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